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Abstract: Using the recently derived NNLO cross sections [1], we provide NNLO+NNLL 
theoretical predictions for top quark pair production based on all the available NNLO 
PDF sets, and compare them with the most precise LHC and Tevatron data. In this 
comparison we study in detail the PDF uncertainty and the scale, mt and dependence 
of the theoretical predictions for each PDF set. Next, we observe that top quark pair 
production provides a powerful direct constraint on the gluon PDF at large x, and include 
Tevatron and LHC top pair data consistently into a global NNLO PDF fit. We then explore 
the phenomenological consequences of the reduced gluon PDF uncertainties, by showing 
how they can improve predictions for Beyond the Standard Model processes at the LHC. 
Finally, we update to full NNLO+NNLL the theoretical predictions for the ratio of top 
quark cross sections between different LHC center of mass energies, as well as the cross 
sections for hypothetical heavy fourth-generation quark production at the LHC. 
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1 Introduction 

The total it production cross section is an important observable at hadron colliders, which 
has been recently computed in full next-to-next-to leading order (NNLO) QCD [1-4]. 
Thanks to a much improved control over higher order terms, including soft gluon emissions 
through next-to-next-to leading log (NNLL) [5-7], scale uncertainties are now controlled 
down to the 2.2 (3) percent level at the Tevatron (LHC), enabling a number of precision 
phenomenology applications to SM and BSM physics. 

In this work we present an in-depth study of the theoretical uncertainties that affect 
the total cross section computed with NNLO+NNLL precision. These uncertainties are 
the parton distributions of the proton, the value of the strong coupling as{Mz), the value 
of the top quark mass mt and the scale uncertainties from missing higher perturbative 
orders. We then compare the theoretical predictions with the most precise available data 
from the Tevatron and the LHC at 7 and 8 TeV. We also provide predictions for LHC 
at 14 TeV, as well as for the ratio of cross sections between 7, 8 and 14 TeV and for 
the production of heavy top-like fermions. Previous phenomenological studies of the total 
tt cross sections, based on different approximations to the full NNLO calculation were 
presented in Refs. [7-14]. 

As a first phenomenological application of the full NNLO calculation we study the im- 
pact of top quark cross section data in the parton distribution analysis. Indeed, top quark 
production is directly sensitive to the large- a; gluon PDF, which at present is affected by 
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substantial uncertainties. In turn, large-x gluons play an important role in theoretical pre- 
dictions of many BSM scenarios like gluino pair production [15], high-mass Kaluza-Klein 
graviton production [16-18], resonances in the invariant mass spectrum [19, 20], quark 
compositeness in inclusive jet and dijet production [21-24] and many others. The availabil- 
ity of the full NNLO calculation makes top quark pair production the only hadron collider 
process that is both sensitive to the gluon and can be consistently included in a NNLO 
PDF fit without any approximations. Hadronic constraints on the gluon PDF are provided 
also by inclusive jet and dijet production [25-28] and isolated photon production [29, 30], 
though these two processes are only known to NLO and affected by substantial scale un- 
certainties.^ 

The focus of this paper is, on the one hand, to provide an up-to-date summary of 
the theoretical uncertainties on the total tt cross section, and on the other hand, to show 
how top quark data can be used to constrain the large-x gluon PDF. Indeed, unlike the 
Tevatron, top quark pair production at the LHC is dominated by gg scattering, thus 
providing a complementary probe of the gluon PDF. As shown in Table 1, at the LHC the 
relative contribution of the gg subprocess is between 85% and 90% depending on the beam 
energy, with qq being about 10-15%, almost the opposite of the Tevatron. 





TeVatron 


LHC 7 TeV 


LHC 8 TeV 


LHC 14 TeV 


99 


15.4% 


84.8% 


86.2% 


90.2% 


qg + Q9 


-1.7% 


-1.6% 


-1.1% 


0.5% 




86.3% 


16.8% 


14.9% 


9.3% 



Table 1. The relative contribution of the various partonic sub-channels to the NNLO+NNLL cross 
section for different colliders and collider energies, computed with the MSTW2008NNLO PDFs. 
We loosely label with qq the sum of all processes without gluons in the initial state. 

To illustrate the range of Bjorken-x's to which the top cross section is sensitive, the 
correlation [32] between the top quark production cross section and the gluon and the up 
quark PDFs is shown in Fig. 1 for the various cases that we will discuss in the paper: 
Tevatron Run II, LHC 7, 8 and 14 TeV. A correlation whose absolute magnitude is close 
to 1 indicates that variations of PDFs with a particular value of x will in turn translate 
into cross-section variations. It is clear from Fig. 1 that for the LHC the top quark cross 
section directly probes the gluon in the range of x between x = 0.1 and x = 0.5, where 
gluon PDF uncertainties are relatively large. 

The outline of this paper is as follows. In Sect. 2 we discuss the settings of the 
calculation and the treatment of the various theoretical uncertainties. In Sect. 3 we provide 
up-to-date predictions for the ti cross section at the Tevatron and LHC and compare with 
the most recent experimental data. In Sect. 4 we quantify the impact of the available 
top data on the gluon PDF, show how it reduces the gluon PDF's large-x uncertainties, 

^Recent progress on the NNLO cross section for jet production was presented in Ref. [31], so in the near 
future it should also be possible to consistently include this process in NNLO PDF fits. 
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Correlation between PDFs and Cross-Section 



Correlation between PDFs and Cross-Section 




Figure 1. The correlation between tire gluon PDF (left plot) and the up quark PDF (right 
plot) from NNPDF2.3 with the total NNLO-fNNLL top quark production cross sections at the 
Tevatron and the LHC for different center of mass energies. The correlations are computed for 
g = mt = 173.3 GeV. 



and study the phenomenological consequences of this improvement. In Sect. 5 we provide 
up-to-date predictions for the ratio of top quark cross sections between 7, 8 and 14 TeV 
and in Sect. 6 we provide predictions for a heavy top-like fermion T. 



2 Settings 

In the following we present the settings of the computation. We use the program Top++, 
v2 . [33] to evaluate the total top quark production cross section at NNLO-I-NNLL accu- 
racy. We use as input the most up-to-date NNLO PDF sets from each PDF collaboration: 
ABMll [34], CTIO [35], HERAPDF1.5 [36], MSTW08 [37] and NNPDF2.3 [38]. The value 
of as{Mz) is set to the preferred value of each group, respectively 0.1134, 0.118, 0.1176, 
0.1171 and 0.118. For NNPDF2.3 we use the set with a maximum of Nf = 5 active fla- 
vors. When providing the predictions for each PDF set, we will use the default in each 
case, while later in Sect. 3.3 we will explore the dependence of the results on the strong 
coupling as ctg is varied. A detailed benchmark comparison of these five NNLO PDF sets 
was recently presented in [39], where the similarities and differences between each of the 
five sets are discussed. 

We consider the following sources of theoretical uncertainties in the top quark pair 
production cross section: 

• Higher perturbative orders. 

The central scales of the NNLO-I-NNLL computation are set to iir = fip = mt- As 
customary, we explore the effect of missing higher perturbative orders by varying the 
scales independently by a factor of two upwards and downwards, with the constraint 
that the ratio of the two scales can never be larger than two. The envelope of the 
resulting cross section defines the scale uncertainty of the computation. As shown in 
Ref. [1], scale variations at the LHC with the full NNLO-I-NNLL result are substan- 
tially smaller that with the NLO or previous approximated NNLO computations. 
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• Parton Distributions. 

We use the corresponding prescription from each group to provide the 68% confidence 
level PDF uncertainties. For the Hessian sets, we use the asymmetric expression 
for PDF uncertainties [40]. The CTIO errors are rescaled by a factor 1.642 since 
they are provided as 90% CL. In the case of HERAPDF1.5, we consider only the 
experimental uncertainties, but not the model and parametrization uncertainties.^ 
Parton distributions are accessed through the LHAPDF interface, version 5.8.9 [44]. 

• Strong coupling constant. 

We assume the 68% CL uncertainty on as{Mz) to be Aa^ = 0.0007, as indicated 
by the 2012 PDG average [41]. For each PDF set, we add in quadrature PDF and 
as uncertainties, except for ABMll, where Og errors are already part of the total 
PDF uncertainty, and thus we avoid double counting. For Hessian PDF sets the 
addition in quadrature of PDF and ag uncertainties is known to be exact [42], while 
for Monte Carlo sets it is a good approximation [43]. The slope of the cross section 
with as for each PDF set is extracted from a linear fit in a common range of as{Mz) 
between 0.116 and 0.120, close to the PDG average and where all PDF groups provide 
predictions. 

• Top quark mass. 

We take as central value mt = 173.3 GeV, with an uncertainty of Arrit = ±1 GeV. 
This is consistent both with the latest Tevatron average,^ mt = 173.20 it 0.87 GeV, as 
well as with the 2012 PDG average [41] of the top quark mass mt = 173.5 it 0.6 it 0.8 
GeV. These masses are to be understood as pole masses. To determine the top-mass 
related uncertainty 6mt of the top cross section, we compute the central value of the 
cross-section for top masses in a 1 GeV range around the central value mt = 173.3. 

Following the recommendations of the Higgs Cross Section Working Group [46] regard- 
ing theoretical uncertainties, we adopt the most conservative method to combine them: we 
add all parametric uncertainties (PDFs, ag and mt) in quadrature and then the total 
parametric error is added linearly to the scale uncertainty to define the total theoretical 
uncertainty. 

3 The top cross-section: measurements, predictions and uncertainties 

With the settings detailed in the previous section, we have computed the cross sections for 
the Tevatron and the LHC and their associated uncertainties. The most precise experi- 
mental measurements of the top cross section are collected in Table 2, while the theoretical 
predictions are collected in Tables 3 to 6. In each case, we provide the NNLO+NNLL 
predictions for the total top quark cross section, for all the PDFs considered, and with the 

^If the HERAPDF1.5 uncertainties for the tt cross section include the model and parametrization un- 
certainties the total PDF uncertainty increases to about 20% [39]. 

^See [45], and the update prepared for the 2013 Rencontres de Moriond QCD, presented by G. Petrillo. 
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Measurement 


o-tf (pb) 


stat. (pb) 


sys. (pb) 


lumi. (pb) 


total (pb) 


Tevatron CDF+DO (Ref. [47]) 


7.65 


± 0.20 


± 0.29 


± 0.22 


7.65 ± 0.42 (5.5%) 


Atlas 7 TeV (Ref. [48]) 


177 


± 3 


+8 
-7 


± 7 


177+10 (+5.6%) 
' -11 (-6.2%) 


CMS 7 TeV (Ref. [49]) 


160.9 


± 2.5 


+5.1 
-5.0 


± 3.6 


160.9 ± 6.6 (4.0%) 


Atlas 8 TeV (Ref. [50]) 


241 


± 2 


± 31 


± 9 


241 ± 32 (13.0%) 


CMS 8 TeV (Refs. [51, 52]) 


227 


± 3 


± 11 


± 10 


227 ± 15 (6.7%) 



Table 2. The experimental measurements from the Tevatron and LHC 7 and 8 TeV used in this 
paper. The statistical, systematic and luminosity uncertainties are added in quadrature for the 
total uncertainty in the last column. See the text for more details. 



various sources of theoretical uncertainties. The default value of as{Mz) from each PDF 
set has been used in the computation of the central predictions. Also included in these 
tables are the best available experimental measurements at the various colliders. 

For the Tevatron, we have used the most up-to-date CDF and DO combinations [47] 
with 8.8 fb-i of data, see Table 2. For LHC 7 TeV, we quote separately the two most 
precise measurements from ATLAS and CIVIS. From ATLAS we use the combination of all 
channels with luminosities between 0.70 and 1.02 fb^^ [48], see Table 2, while for CMS 
we use the single most precise measurement^, obtained in the dilepton channel with a 
luminosity of 2.03 fb"^ [49], also given in Table 2. Note also that a preliminary combination 
of all channels between ATLAS and CIViS for an integrated luminosity of 1.0 fb~^ has been 
presented in [53], yielding a{tt) = 173.3 ± 10.1 pb. For LHC 8 TeV, there are three 
available measurements, two from CIVIS, in the di-lepton channel, with 2.4 fb~^ [51] and 
in the lepton+jets channel with 2.8 fb~^ [52], and one from ATLAS in the lepton+jets 
channel with, 5.8 fb^-"^ [50]. We use the CMS dilepton and the ATLAS lepton+jets results, 
collected in Table 2. 

In Fig. 2 we show the predictions from each PDF set compared to experimental data. 
The X-axis is the reference value of for each PDF set. The inner error bars correspond 
to the linear sum of scale and PDF uncertainties, while the outer error bar is the total 
theoretical uncertainty including also the contributions from 6mt and Sag- 

In the remainder of this section we comment on various features of the above results. 

3.1 Comparison betv^een the various PDF sets 

We notice that the results for NNPDF2.3, CTIO, MSTW08 and HERAPDF1.5 are all close 
to each other, both in central values and theoretical uncertainties. On the other hand, the 
ABlVIll predictions are lower compared to other sets and to experimental data. Fig. 1 
suggests that similarities or differences in the total cross section can be understood in 
terms of the large-x behavior of the gluon PDF, together with the value of as used by each 
PDF set. 

*In this case, given that the dependence of the measured total cross section with rrit is provided, we 
have rescaled the central value of the data to the same value of the top quark mass used for the theory 
computations, nit ~ 173.3 CeV. 
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Tevatron Run II 



PDF set 


CTtt (pb) 


f^scalc 


(pb) 




(pb) 




(pb) 




(pb) 


(5tot (pb) 


ABMll 


6.869 


+ 0.104 
-0.174 


( + 1.5%) 
(-2.5%) 


+0.157 
-0.157 


(+2.3% ) 
(-2.3%) 


+0.000 
-0.000 


(+0.0%) 
(-0.0%) 


+0.207 
-0.201 


(+3.0%) 
(-2.9%) 


+0.364 (+5.3%) 
-0.429 (-6.2%) 


CTIO 


7.395 


+0.116 
-0.210 


( + 1.6%) 
(-2.8%) 


+0.270 
-0.203 


(+3.6%) 
(-2.7%) 


+0.136 
-0.136 


(+1.8%) 
(-1.8%) 


+0.235 
-0.227 


(+3.2%) 
(-3.1%) 


+0.499 (+6.7%) 
-0.544 (-7.4%) 


HER A 1.5 


7.624 


+0.116 
-0.074 


(+1.5%) 
(-1.0%) 


+0.134 
-0.154 


(+1.8%) 
(-2.0%) 


+0.098 
-0.098 


(+1.3%) 
(-1.3%) 


+0.241 
-0.233 


(+3.2%) 
(-3.1%) 


+0.409 (+5.4%) 
-0.370 (-4.9%) 


MSTW08 


7.164 


+0.110 
-0.200 


( + 1.5%) 
(-2.8%) 


+0.169 
-0.122 


(+2.4%) 
(-1.7%) 


+0.088 
-0.088 


(+1.2%) 
(-1.2%) 


+0.228 
-0.220 


(+3.2%) 
(-3.1%) 


+0.391 (+5.5%) 
-0.475 (-6.6%) 


NNPDF2.3 


7.258 


+0.117 
-0.202 


( + 1.6%) 
(-2.8%) 


+0.121 
-0.121 


(+1.7%) 

(-1.7%) 


+0.090 
-0.090 


(+1.2%) 
(-1.2%) 


+0.229 
-0.221 


(+3.1%) 
(-3.0%) 


+0.390 (+5.4%) 
-0.469 (-6.5%) 


CDF+DO 


7.65 










± 0.42 ( 5.5%) 



Table 3. The NNLO+NNLL predictions for the total top quark pair cross-section at the Tevatron 
Run II, for ah the PDFs considered, and with the various sources of theoretical uncertainties. The 
default value of as{Mz) from each collaboration has been used in the computation of the central 
predictions. The four theoretical uncertainties are combined into a total theory error as discussed 
in the text. The lower row shows the best available experimental measurement. 



The agreement between CT, HERAPDF, MSTW and NNPDF can be traced back to 
(a) a similar default value of Ug used and (b) a similar large-x gluon PDF [39]. On the 
other hand, regarding the differences between ABMll and the other sets, we note (a) the 
smaller value of ag used by ABMll and (b) the softer large-x gluon PDF in the region 
relevant for top quark production [39]. As shown below in Fig. 3, using a value for Us closer 
to the PDG average would improve the agreement of the ABMll set both with other PDF 
sets and with experimental data, though the ABMll predictions at the LHC are still lower 
than that of other sets even for a common value. 

To understand the differences in the large-x gluon in ABMll as compared to other 
sets, Refs. [57, 58] have suggested that the reason is that, keeping everything else fixed, 
the use of a fixed- flavor number scheme in a fit to deep- inelastic data (such as ABMll) 
leads to a softer large-x gluon as compared to the gluon obtained in PDF fits based on 
variable-flavor number schemes [54-56]. On the other hand, the differences between the 
various implementations of variable-flavor number schemes translate into much reduced 
differences in the PDFs [57] and thus into the top cross sections. 

In order to quantify the agreement of data and theory, we compute the for the 
central values of the five PDF sets under consideration and for the five available data 
points of Table 2: one for the Tevatron, two for LHC 7 TeV and two for LHC 8 TeV. For 
each of the five PDF sets studied, the is defined as 



\ tt tt 



1=1 (>tot 

with A'dat = 5. We consider all five data points as fully uncorrelated.^ In addition. 



^This is a necessary approximation, but as discussed in [53] there is some degree of correlation between 
the ATLAS and CMS data, as there is between LHC data at different cm. energies. 
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LHC 7 TeV 



PDF set 


CTtt (pb) 


(Jscalo (pb) 


5pDF (pb) 




(pb) 




(pb) 


(5tot (pb) 


ABMll 


135.8 


+3.5 
-4.2 


(+2.6%) 
(-3.1%) 


+6.4 
-6.4 


(+4.7%) 
(-4.7%) 


+0.0 
-0.0 


( + 0.0%) 

(-0.0%) 


+4.3 
-4.2 


(+3.2%o) 
(-3.1%) 


+ 11.2 (+8.270) 
-11.8 (-8.7%) 


CTIO 


172.5 


+4.6 
-6.0 


(+2.7%) 
(-3.5%) 


+8.0 
-6.5 


(+4.6%) 
(-3.8%) 


+3.7 
-3.7 


(+2.2%) 
(-2.2%) 


+5.3 
-5.1 


(+3.1%) 
(-3.0%) 


+ 14.9 (+8.6%) 
-15.0 (-8.7%) 


HER A 1.5 


177.2 


+4.8 
-4.2 


(+2.7%) 
(-2.3%) 


+4.0 
-6.4 


(+2.3%) 
(-3.6%) 


+3.0 
-3.0 


(+1.7%) 
(-1.7%) 


+5.4 
-5.2 


(+3.1%) 
(-2.9%) 


+ 12.2 (+6.9%) 
-12.9 (-7.3%) 


MSTW08 


172.0 


+4.4 
-5.8 


(+2.6%) 
(-3.4%) 


+4.7 
-4.7 


(+2.7%) 
(-2.7%) 


+2.9 
-2.9 


(+1.7%) 
(-1.7%) 


+5.3 
-5.1 


(+3.1%) 
(-3.0%) 


+ 12.1 (+7.0%) 
-13.4 (-7.8%) 


NNPDF2.3 


172.7 


+4.6 
-6.0 


(+2.7%) 
(-3.5%) 


+5.2 
-5.2 


(+3.0%) 
(-3.0%) 


+2.7 
-2.7 


(+1.6%) 

(-1.6%) 


+5.3 
-5.2 


(+3.1%) 
(-3.0%) 


+ 12.5 (+7.2%) 
-13.7 (-8.0%) 


ATLAS 
CMS 


177 
160.9 










+10 (+5.6%) 
-11 (-6.2%) 

± 6.6 ( 4.0%) 



Table 4. Same as Table 3 for LHC 7 TeV. 



LHC 8 TeV 



PDF set 


Ttt (pb) 




(pb) 


5pDF (pb) 




(pb) 




(pb) 


Stot 


(pb) 


ABMll 


198.6 


+5.0 
-6.2 


+2.5%) 
-3.1%) 


+8.5 
-8.5 


(+4.3%) 
(-4.3%) 


+0.0 
-0.0 


(+0.0%) 
(-0.0%) 


+6.1 
-5.9 


(+3.1%) 
(-3.0%) 


+ 15.5 
-16.6 


( + 7.8%) 
(-8.4%) 


CTIO 


246.3 


+6.4 
-8.6 


+2.6%) 
-3.5%) 


+ 10.1 
-8.2 


(+4.1%) 
(-3.3%) 


+4.9 
-4.9 


(+2.0%) 
(-2.0%) 


+7.4 
-7.1 


(+3.0%) 
(-2.9%) 


+ 19.8 

-20.5 


(+8.1%) 
(-8.3%) 


HERA1.5 


252.7 


+6.5 
-5.9 


+2.6%) 
-2.3%) 


+5.4 
-8.6 


(+2.1%) 
(-3.4%) 


+4.0 
-4.0 


(+1.6%) 
(-1.6%) 


+7.5 
-7.3 


(+3.0%) 
(-2.9%) 


+ 16.6 
-17.8 


(+6.6%) 
(-7.1%) 


MSTW08 


245.8 


+6.2 
-8.4 


+2.5%) 
-3.4%) 


+6.2 
-6.2 


(+2.5%) 
(-2.5%) 


+4.0 
-4.0 


(+1.6%) 
(-1.6%) 


+7.4 
-7.1 


(+3.0%) 
(-2.9%) 


+ 16.6 

-18.7 


(+6.8%) 
(-7.6%) 


NNPDF2.3 


248.1 


+6.4 
-8.7 


+2.6%) 
-3.5%) 


+6.6 
-6.6 


(+2.7%) 
(-2.7%) 


+3.7 
-3.7 


(+1.5%) 
(-1.5%) 


+7.5 
-7.2 


(+3.0%) 
(-2.9%) 


+ 17.1 
-19.1 


(+6.9%) 
(-7.7%) 


ATLAS 


241.0 


















± 32.0 


( 13.3%) 


CMS 


227.0 


















± 15.0 ( 6.6%) 



Table 5. Same as Table 3 for LHC 8 TeV. 



LHC 14 TeV 



PDF set 


att (pb) 




(pb) 


<5pDF 


(pb) 




(pb) 




(pb) 


(5tot 


(pb) 


ABMll 


832.0 


+ 18.7 
-27.4 


+2.2%) 
-3.3%) 


+25.1 
-25.1 


(+3.0%) 
(-3.0%) 


+0.0 
-0.0 


(+0.0%) 
(-0.0%) 


+23.3 
-22.5 


(+2.8%) 
(-2.7%) 


+52.9 
-61.1 


(+6.4%) 
(-7.3%) 


CTIO 


952.8 


+23.3 
-34.5 


+2.4%) 
-3.6%) 


+22.4 
-19.9 


(+2.3%) 
(-2.1%) 


+ 14.0 

-14.0 


(+1.5%) 

(-1.5%) 


+26.1 
-25.2 


(+2.7%) 
(-2.6%) 


+60.3 
-69.5 


(+6.3%) 
(-7.3%) 


HER A 1.5 


970.5 


+22.1 
-22.0 


+2.3%) 
-2.3%) 


+ 15.3 
-25.7 


(+1.6%) 
(-2.6%) 


+ 12.8 
-12.8 


(+1.3%) 
(-1.3%) 


+26.4 
-25.6 


(+2.7%) 
(-2.6%) 


+55.2 
-60.5 


(+5.7%) 
(-6.2%) 


MSTW08 


953.6 


+22.7 
-33.9 


+2.4%) 
-3.6%) 


+16.2 
-17.8 


(+1.7%) 
(-1.9%) 


+ 12.8 
-12.8 


(+1.3%) 
(-1.3%) 


+26.1 
-25.3 


(+2.7%) 
(-2.7%) 


+56.3 
-66.8 


(+5.9%) 
(-7.0%) 


NNPDF2.3 


977.5 


+23.6 
-35.4 


+2.4%) 
-3.6%) 


+ 16.4 
-16.4 


(+1.7%) 
(-1.7%) 


+ 12.2 
-12.2 


(+1.3%) 

(-1.3%) 


+26.9 
-26.1 


(+2.8%) 
(-2.7%) 


+57.4 
-68.5 


(+5.9%) 
(-7.0%) 



Table 6. Same as Table 3 for LHC 14 TeV. 



to quantify the overall consistency of theoretical predictions with experimental data, we 



TeVatron 



LHC 7 TeV 



# ABM1 1 

■ CT10 

A HERAPDF 

•k MSTW2008 

y NNPDF2.3 



Tevatron Run-ll 



0.112 0.113 0.114 0.115 0.116 0.117 0.118 0.119 0.12 
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Figure 2. The best predictions from eaclr PDF set compared to experimental data, as a function 
of the default as{Mz) value. The inner error bar corresponds to the linear sum of PDF and scale 
uncertainties, while the outer error bar is the total theoretical uncertainty, computed as described 
in the text. 



define a pull estimator as follows, 



(cxp) _ (th)^ 2 
"tt "tt 



N^t^ .(cxp)2 .(th)2 ' (3-2) 



^tot I "^tot 

where now in the denominator we add in quadrature the experimental and theoretical total 
uncertainties. For simplicity we have symmetrized the total theoretical systematic error 
for each PDF set. 

The results for both the ^i^d the pull are summarized in Table 7. We provide both 
the total the individual contribution from the data at the Tevatron, LHC 7 and 

8 TeV. As we can see, most PDF sets provide a good description of the total top quark 
cross section data, with x^/^dat ~ 1- The pull also shows the good consistency between 
experimental data and theory prediction for most of the PDF sets considered. 

3.2 Uncertainty due to the value of mt 

The uncertainty on the value of the top quark mass is now a substantial fraction of the 
overall systematic uncertainty in the total tt cross section. In Table 8 we compare, for the 
NNPDF2.3 set, the total theory uncertainty with and without including the uncertainty 
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Xtcv 




xiiics 


Xtot 


X?ot/^dat 


P 


AMBll 


3.5 


31.4 


5.3 


40.2 


8.0 


3.2 


CTIO 


0.4 


3.3 


1.7 


5.3 


1.1 


0.3 


HERAPDF15 


0.0 


6.1 


3.1 


9.2 


1.8 


0.5 


MSTW08 


1.3 


3.1 


1.6 


6.0 


1.2 


0.4 


NNPDF2.3 


0.9 


3.4 


2.0 


6.3 


1.3 


0.4 



Table 7. The between data and theory, Eq. (3.1), for all the five PDF sets, both for the total 
dataset and spht into colliders. Let us recall that the standard deviation of the distribution 
for A^dat = 5 data points is \/2Ndat = 3.1. The last column show the pull Eq. (3.2) between 
experimental data and theory predictions. 

due to the top quark mass. We see that at the LHC, an uncertainty of 5mt = ±1 GeV 
translates into an increase between 1 and 1.5% of the total theory uncertainty, and into a 
somewhat larger increase at the Tevatron. Given that PDFs, scale and mt uncertainties 
are all now of similar size, the total theory error would be only slightly reduced if one 
assumed that the uncertainty due to mt could be completely neglected, for example after 
more precise measurements of this parameter at the LHC. 



Collider 


o-tt (pb) 


^PDF+scales+Qs (pb) 


<5tot (pb) 


Tevatron 
LHC 7 TcV 
LHC 8 TeV 
LHC 14 TeV 


7.258 
172.7 
248.1 
977.5 


+0.267 (+3.7%) 
-0.3.52 (-4.9%) 

+ 10.4 (+6.0%) 
-11.8 (-6.8%) 

+ 14.0 (+5.6%) 
-16.2 (-6.5%) 

+44.1 (+4.5%) 
-55.8 (-5.7%) 


+0.390 (+5.4%) 
-0.469 (-6.5%) 

+ 12.5 (+7.2%) 
-13.7 (-8.0%) 

+ 17.1 (+6.9%) 
-19.1 (-7.7%) 

+57.4 (+5.9%) 
-68.5 (-7.0%) 



Table 8. The NNLO+NNLL predictions for the total top quark pair cross-section at the Tevatron 
and the LHC, with NNPDF2.3 as input. The third column shows the theoretical uncertainty when 
5mt is not taken into account, while the last column is the total theoretical uncertainty (the same 
as in Tables 3 to 6). At the LHC, an uncertainty of 5rat ~ ±1 GcV translates into an increase 
between 1 and 1.5% of the total theory uncertainty, and a bit more at the Tevatron. 

3.3 Uncertainty due to the value of Ug 

In Fig. 3 we show the dependence with as{Mz) of the NNLO+NNLL cross sections for 
each of the various PDF sets. We use consistently the same value of as{Mz) in the partonic 
cross sections and in the PDFs, using all the range of varying PDFs provided by each 
group. We also include for reference the best experimental data. The slope with as is 
similar for each of the PDF sets. One can see that a variation of Aq^ by 0.001 increases 
the cross section by about 0.13 pb at the Tevatron, 4 pb at LHC 7 TeV, 6 pb at 8 TeV and 
20 pb at 14 TeV, with small differences between PDF sets. 

We notice that the Ug dependence of the total cross section reported in our tables 
is slightly larger than what one would naively estimate based on the power counting of 
the partonic cross section. The reason is that, in the range of x relevant for top quark 
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Figure 3. The theoretical predictions for tlie various PDF sets studied in this work as a function 
of the strong coupHng constant as{Mz), compared to the best available experimental data. From 
top to bottom we show Tevatron, LHC 7 TeV, 8 TeV and 14 TeV. 



production, the value of the gluon density is larger for PDF fits with a larger as- This is 
shown in Fig. 4, for the specific case of NNPDF2.3, but the behavior is similar for other 
sets: larger ag leads to smaller g[x) in the small/medium-x region that is controlled by 
deep-inelastic HERA data, while DGLAP evolution and the momentum sum rule balance 
this reduction with an increase at larger x. 

The net effect of this positive correlation is an enhanced sensitivity to of the ti 
cross section, about 30% larger than what one would have obtained by keeping the PDFs 
fixed while changing a^. This fact, together with the small, 2-3%, scale uncertainty of 
the full NNLO+NNLL result, suggest that the top production cross section could provide 
a useful independent determination of as{Mz)- Such extraction would be analogous to 
the determinations of rrit from the same total cross section [12-14, 59-61]. The strong 
coupling can also be determined at the Tevatron and the LHC from jet production [62-64], 
but results are affected by sizable scale uncertainties from the lack of a complete NNLO 
calculation as well as by non-perturbative corrections. A first determination of ag from 
the ti cross section, using previous approximate NNLO predictions, has been performed by 
the CMS Collaboration [65] . 
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Ratio to NNPDF2.3 NNLO (ag = 0.119), = mf 




Figure 4. The gluon density in NNPDF2.3 fits obtained for different input values of as{Mz), 
evaluated at = mf. Results are shown as the ratio with respect the gluon PDF obtained for 
a,(Mz) =0.119. 



4 Top quark data constraints on the gluon PDF 

The gluon PDF is one of the worse known partonic distributions. Deep-inelastic data 
constrain the gluon only indirectly, via scaling violations, and direct information comes only 
from collider inclusive jet production data [25-28]. Recently, the use of isolated photon data 
has also been advocated to pin down the gluon [29]. However, direct constrains at hadron 
collider from photon and jet data are affected by substantial scale uncertainties due to the 
missing full NNLO result, and are complicated by various non-perturbative uncertainties. 
The availability of a full NNLO calculation makes the total top pair cross section the only 
collider observable which, at present, is both directly sensitive to the gluon PDF and can be 
consistently included in a NNLO QCD analysis. The fact that non-perturbative corrections 
are much reduced in the total top quark production cross-section as compared to photons 
and jets is another good motivation to use tops as probes of the gluon PDF. The possibility 
of using top cross section data in PDF analysis was already discussed at the qualitative 
level in Ref. [66]. More recently, Ref. [13] provided a first estimate of the impact of top 
quark data on the gluon PDF based on approximate NNLO results. 

In Sect. 3 we have shown at the qualitative level that available data already discrimi- 
nates between different PDF sets, that is, between different large- a; gluon PDFs. Now we 
will be more quantitative and determine if available data can help in reducing the gluon 
uncertainties within a single PDF set. In order to quantify this impact, we will use the 
Bayesian PDF reweighting method of Refs. [67, 68] on the NNPDF2.3 set. The same study 
could be carried out with Hessian PDF sets as discussed in [69]. In principle, one could also 
use the top quark differential distributions data from ATLAS and CMS [70, 71]. However, 
these data are less precise than the total cross section measurement and the corresponding 
theoretical predictions are currently only available at NLO.® 

®For some specific differential distributions, results at NLO supplemented by threshold resummation are 
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Figure 5. Left plot: the ratio of the NNPDF2.3 NNLO gluon PDF at = iqO GeV^ between 
the defauh fit and after including the Tcvatron and LHC top quark cross section data. Right plot: 
the relative reduction of PDF uncertainties thanks to the inclusion of top data in the PDF fit. 



Therefore, we have included the N^at = 5 experimental data points available from the 
Tevatron and the LHC into the NNPDF2.3 NNLO fit. The definition of that we use 
is Eq. (3.1). The effective number of replicas after reweighting (the exponential of the 
Shannon entropy) is Ncg = 86, out of the starting sample of 100 replicas, indicating the 
moderate constraining power of the data. 

The results of adding the top quark data into the gluon PDF are shown in Fig. 5. We 
show the NNPDF2.3 NNLO gluon at = 100 GeV^, in the default fit and after including 
the Tevatron and LHC top quark production cross section data. We observe that the large- 
X gluon PDF uncertainties decrease. This is expected since in that region the correlation 
between the gluon and the top-quark cross section is maximal (see Fig. 1). In Fig. 5 we 
also show the relative reduction in PDF uncertainties from the addition of top data. This 
reduction coincides, within the finite statistics of the original NNPDF Monte Carlo sample, 
with the correlation profile of Fig. 1, and in particular confirms that the top quark data 
have a small impact below x ~ 0.1 or so, as well as for very large values of x. Therefore, 
we conclude that the available top data can already help reduce the uncertainties on the 
gluon PDF by a factor of up to 20%, and are thus an important ingredient to future global 
PDF analyses. On the other hand, the impact of top data on the quark PDFs is essentially 
negligible. 

It is interesting to study the modifications of the theory predictions after the top 
quark data have been added into the NNPDF2.3 fit. In Table 9 we show the tt cross 
section for NNPDF2.3, comparing the default prediction with the predictions after adding 
different subsets of the top quark data. We show only the entries which correspond to pure 
predictions. By including top data from lower energy colliders, we can provide arguably 
the most accurate theoretical prediction for the total ti cross section at higher energies, 
given that PDF uncertainties will be reduced in the same kinematical range from lower 
energy data. 

These predictions are collected in Table 9. As an illustration, the NNPDF2.3 prediction 

available, see [9, 72] and references therein. 
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including Tevatron and LHC 7 top data would be the best available theory prediction for 
LHC 8 TeV. Note that not only PDF uncertainties are reduced, but that also the central 
value is shifted to improve the agreement with the experimental data. As can be seen, the 
precise 7 TeV data carry most of the constraining power, though of course improved power 
of the 8 TeV data will be provided with the analysis of the full 2012 dataset. 



Collider 


Ref 


Ref+TcV 


Ref +TeV+LHC7 


Ref+TcV+LHC7+8 


Tevatron 


7.26 ± 0.12 








LHC 7 TeV 
LHC 8 TeV 
LHC 14 TeV 


172.5 ± 5.2 
247.8 ± 6.6 
976.5 ± 16.4 


172.7 ± 5.1 
248.0 ± 6.5 
976.2 ± 16.3 


245.0 ± 4.6 
969.8 ± 12.0 


969.6 ± 11.6 



Table 9. The tt cross section CTj^, in picobarns, for the NNPDF2.3 NNLO set, together with the 
associated PDF uncertainties. We show both the reference predictions, and the predictions once 
Tevatron, LHC7 and LHC8 data are added sequentially to the fitted dataset. We only show the 
predictions corresponding to beam energies whose data have not been included in the fit. The 
fourth column corresponds to the best theory prediction for LHC 8 TeV, while the last column is 
the best theory prediction for LHC 14 TeV. 

Then in Table 10 we provide NNPDF2.3 compared to the top quark data, before 
adding any data, after adding all Tevatron and LHC data and adding only the Tevatron 
and LHC 7 TeV data points. The slight improvement of an already good quantitative 
description can be seen. As expected, the agreement of the prediction with LHC8 data, 
when only Tevatron and LHC7 data are used, is a non-trivial consistency check of the 
whole procedure. 



Collider 


(Total, iVdat = 5) 


(LHC 8 TeV, Vdat = 2) 


NNPDF2.3 


6.28 


1.64 


NNPDF2.3 + TeV,LHC data 


4.88 


1.24 


NNPDF2.3 + TeV,LHC7 data 


4.87 


1.24 



Table 10. The NNPDF2.3 compared to the hadron collider top quark production data, before 
adding any data (first line), after adding all Tevatron and LHC data (second line) and adding only 
the Tevatron and LHC 7 TeV data points (third line). Note that the is not normalized by the 
number of data points. 



Given that the constraints from top quark data in a global PDF fit such as NNPDF2.3 
are already substantial, we expect even larger constrains in PDF sets based on reduced 
datasets. To quantify the impact of top data into a DIS~only PDF fit, we have performed 
a similar analysis as the one with NNPDF2.3, but now starting from a PDF set based on a 
reduced dataset, the NNPDF2.1 HERA-only set [73, 74]. As the name indicates, this PDF 
set includes only HERA data, and thus is affected by larger PDF uncertainties than the 
NNPDF global fits, in particular for the large-x gluon. 

In Fig. 6 we compare first of all the gluon from the NNPDF2.1 HERA-only fit with 
that of the NNPDF2.1 global fit, to show the large differences in PDF uncertainties due to 
the reduced dataset in the former case. Then we show the improvements in the gluon PDF 
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in the HERA-only fit after the addition of the top quark data. It is clear that the impact 
is substantial. Remarkably, the top quark data bring the gluon from the HERA-only fit 
closer to the gluon from the global fit. 



Ratio to NNPDF2.1 NNLO HERA-only, a^.O.lig Ratio to NNPDF2.1 NNLO HERA-only, aj, = 0.118 




Figure 6. Left plot: the NNPDF2.1 HERA-only gluon PDF, at = iqo GeV^, compared to the 
reference global NNPDF2.1 gluon. Right plot: the NNPDF2.1 HERA-only gluon before and after 
including the LHC top quark production cross sections. 

4.1 Impact on predictions for BSM particle production 

Many scenarios of BSM physics predict the production of massive final states in gluon- 
initiated processes. The improvement in the large-x gluon PDF uncertainties seen in Fig. 5 
implies a similar improvement in any cross section that is dominated by initial-state gluons 
in a similar kinematical region. Production kinematics determines that any high-mass 
final state that is gluon initiated will benefit from the improvement in PDF uncertainties. 
Therefore, now we explore in two cases the phenomenological implications of the improved 
large-x gluon: the production of heavy Kaluza-Klein resonances in warped extra dimensions 
scenarios, and the PDF uncertainties in the high invariant mass distribution of top quark 
pairs. 

To begin with, we have considered Kaluza-Klein massive graviton production in warped 
extra dimension scenarios [17], in particular in the so-called bulk models in which the cou- 
pling of the graviton to fermions is suppressed, and thus production is driven by gluon-gluon 
annihilation [18] . We have computed the cross sections for Randal-Sundrum graviton pro- 
duction at the LHC 8 TeV at leading order using the MadGraphS program [75] , for a range 
of values of the graviton mass Mq- The improvement in theoretical (PDF) uncertainties 
in the graviton production cross section thanks to top quark data is shown in Fig. 7. We 
see that, to begin with, PDF uncertainties are large, almost 40% for Mq = 3 TeV, and 
growing as we approach the kinematic threshold, reflecting our lack of knowledge of the 
large-x gluon PDF. The top quark data reduce moderately the production uncertainties. Of 
course, future, more precise top data will render these constraints more stringent, which in 
turn translate into the possibility of better characterizing eventual BSM high mass particles 
that the LHC could find. Similar improvements are found for 14 TeV. 
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Ratio to NNPDF2.3 NNLO 




Figure 7. The PDF uncertainties for tlie production of Randall-Sundrum Kaluza-Klcin gravitons 
at tlic LHC 8 TeV, with NNPDF2.3 before and after including tlic top quark data, as a function 
of the graviton mass Mq. We have assumed that the graviton couples only to gluons. The cross 
section has been computed at leading order with MadGraphS. 



The second example is the high mass tail of the invariant mass distribution in top- 
quark production, relevant for many BSM searches [19, 20], which is also related with 
the top quark forward-backward asymmetry measured at the Tevatron. One example is 
the searches for heavy resonances that decay into a top-antitop pair [76-79]. To study 
the impact of the improvement of gluon PDF uncertainties there, we have computed top 
quark pair production with the aMCSNLO program [80-82] , at NLO matched to the Herwig6 
shower [83]. The renormalization and factorization scales are set equal to the sum of 
transverse masses of all final-state particles. Using NNPDF2.3, we have evaluated the 
cross section with a cut in the minimum invariant mass of the ti pair, Mtt, before and 
after including the top quark data into the PDF fit. We show the results in Fig. 8, for the 
absolute cross section (left) and for the relative scale and PDF uncertainties (right) as a 
function of the minimum value of Mu allowed. 

As we can see, the addition of the total ti cross section data reduces the PDF uncer- 
tainty in the high-end tail of the ti mass distribution. It must be remarked that the data in 
this high- mass tail are only a negligible fraction of the total ti cross section, and therefore 
do not play any role in the PDF fit itself. The PDF fit including a{ti) reduces the gluon 
uncertainty at large x only because of the overall constraints on the PDF evolution, which 
correlate the x behavior in the x ~ 0.1 region (which dominates the total production cross 
section) and the large- a; region, which is relevant to the high-mass behavior. 

As in the case of dijet cross sections, we expect that rate measurements in kinematical 
regions where, for example, the ti system has a large rapidity, can be used to further improve 
the knowledge of large-x gluons, and improve even more the precision of predictions for 
the production of large-mass objects in (/^-initiated channels. 
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Figure 8. Left plot: the tt cross section above a minimum value of the tt invariant mass Mu, 
computed with aMCatNLO with the Herwig parton shower, and the NNPDF2.3 as input set, at the 
LHC 8 TeV. The error band is the PDF uncertainty band. Right plot: the relative scale and PDF 
uncertainties for the cross section of top quark pairs with invariant mass above Mu.min, at LHC 
8 TeV. In both cases we show the predictions with NNPDF2.3 before and after including the top 
quark data into the PDF fit. 



5 Cross-section ratios between different LHC beam energies 

The measurement of cross-section ratios between different center of mass energies at the 
LHC has two main motivations [84]: first, they are interesting for precision SM studies and 
second, they have the potential to enhance the possible BSM sensitivity of absolute cross 
sections. In Ref. [84], results for top quark cross sections based on NLO+NNLL theory were 
provided; here we update them to NNLO+NNLL and compare them with experimental 
data. 

First of all we show the correlations between PDFs and the cross-section ratios in Fig. 9. 
We see that the PDF correlation is approximately the inverse of that of the absolute cross 
sections, shown in Fig. 1. The reason for this anti-correlation at large x is that when going 
to higher energies, the average probed values of x are smaller. 



Correlation between PDFs and Cross-Section 



Correlation between PDFs and Cross-Section 
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Figure 9. The correlation between the gluon PDF (left plot) and the up quark PDF (right plot) 
with the top quark production cross-section ratios between different LHC center of mass energies. 
The correlations are computed for Q = mt = 173.3 GeV. 
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Using the same settings as before, we now provide NNLO+NNLL predictions for the 
cross-section ratios in Tables 11 and 12. We note that all systematic theoretical uncer- 
tainties are small: in the ratio between 8 and 7 TeV the total theory uncertainty is at the 
permille level, while in the ratio between 14 and 8 TeV, the total error on the ratio is at 
most around 2%. This is dominated by PDF uncertainty, which is four times larger than 
the combined scale, ag and rrit uncertainties. Notice the great improvement due to the 
NNLO result: in the NLO+NNLL analysis of top cross-section ratios of Ref. [84], scale 
and PDF uncertainties were of similar size. We note that, as can be seen from Table 12, 
in the case of the 14 over 8 TeV ratios the differences between PDF sets are at the level of 
10%. Therefore, such measurement would provide powerful PDF constraints even within 
the foreseeable experimental accuracies. 



LHC 8 TeV / 7 TeV ratio 



PDF set 




f^scalc 


<5pDF 




(Pb) 


Suit 


Stot 


ABMll 


1.463 


+0.001 
-0.002 


(+0.1%) 
(-0.1%) 


+0.00C 
-0.006 


(+0.4%) 
(-0.4%) 


+0.000 
-0.000 


(+0.0%) 
(-0.0%) 


+0.001 
-0.001 


(+0.1%) 
(-0.1%) 


+0.007 (+0.5%) 
-0.008 (-0.5%) 


CTIO 


1.428 


+0.001 

-0.001 


(+0.1%) 
(-0.1%) 


+0.008 
-0.010 


(+0.5%) 
(-0.7%) 


+0.002 
-0.002 


(+0.2%) 
(-0.2%) 


+0.001 
-0.001 


(+0.1%) 
(-0.1%) 


+0.009 (+0.6%) 
-0.011 (-0.8%) 


HER A 1.5 


1.426 


+0.001 
-0.002 


(+0.0%) 
(-0.1%) 


+0.003 
-0.003 


(+0.2%) 
(-0.2%) 


+0.001 
-0.001 


(+0.1%) 
(-0.1%) 


+0.001 
-0.001 


(+0.1%) 
(-0.1%) 


+0.004 (+0.3%) 
-0.005 (-0.4%) 


MSTW08 


1.429 


+0.001 
-0.001 


(+0.1%) 
(-0.1%) 


+0.004 
-0.004 


(+0.2%) 
(-0.2%) 


+0.001 
-0.001 


(+0.1%) 
(-0.1%) 


+0.001 
-0.001 


(+0.1%) 
(-0.1%) 


+0.005 (+0.3%) 
-0.005 (-0.3%) 


NNPDF2.3 


1.437 


+0.001 
-0.001 


(+0.1%) 
(-0.1%) 


+0.006 
-0.006 


(+0.4%) 
(-0.4%) 


+0.001 
-0.001 


(+0.1%) 
(-0.1%) 


+0.001 
-0.001 


(+0.1%) 
(-0.1%) 


+0.007 (+0.5%) 
-0.007 (-0.5%) 


ATLAS 


1.36 


















± 0.11 ( 8%) 


CMS 


1.40 


















± 0.08 ( 6%) 



Table 11. The NNLO+NNLL predictions for the ratio of top quark cross section at the LHC 
between 8 and 7 TeV, with all the PDFs considered, and with the various sources of theoretical un- 
certainties. The default value of as{Mz) from each collaboration has been used in the computation 
of the central predictions. The lower row shows our estimate of the best available experimental mea- 
surement. The total theoretical uncertainty is the linear sum of scale and parametric uncertainties, 
as discussed in the text. 



LHC 14 TeV / 8 TeV ratio 



PDF set 


CTtt 


f^scalc 


<5pDF 


5a, 


(pb) 




Stot 


ABMll 


4.189 


+0.008 
-0.016 


(+0.2%) 
(-0.4%) 


+0.057 
-0.057 


(+1.4%) 
(-1.4%) 


+0.000 (+0.0%) 
-0.000 (-0.0%) 


+0.012 
-0.012 


(+0.3%) 
(-0.3%) 


+0.067 (4 
-0.074 (- 


-1.6%) 
-1.8%) 


CTIO 


3.869 


+0.006 
-0.009 


(+0.2%) 
(-0.2%) 


+0.068 
-0.088 


(+1.8%) 
(-2.3%) 


+0.020 
-0.020 


(+0.5%) 
(-0.5%) 


+0.010 
-0.010 


(+0.2%) 
(-0.2%) 


+0.077 (4 
-0.100 (- 


-2.0%) 
-2.6%) 


HER A 1.5 


3.841 


+0.005 
-0.012 


(+0.1%) 
(-0.3%) 


+0.033 
-0.025 


(+0.9%) 
(-0.7%) 


+0.010 
-0.010 


(+0.3%) 
(-0.3%) 


+0.009 
-0.010 


(+0.2%) 
(-0.2%) 


+0.041 (4 
-0.041 (- 


-1.1%) 
-1.1%) 


MSTW08 


3.880 


+0.006 
-0.009 


(+0.2%) 
(-0.2%) 


+0.036 
-0.036 


(+0.9%) 
(-0.9%) 


+0.011 
-0.011 


(+0.3%) 
(-0.3%) 


+0.010 
-0.010 


(+0.2%) 
(-0.2%) 


+0.045 (4 
-0.048 (- 


-1.2%) 
-1.2%) 


NNPDF2.3 


3.940 


+0.006 
-0.010 


(+0.2%) 

(-0.3%) 


+0.048 
-0.048 


(+1.2%) 
(-1.2%) 


+0.009 
-0.009 


(+0.2%) 
(-0.2%) 


+0.010 
-0.010 


(+0.3%) 
(-0.3%) 


+0.056 (4 
-0.060 (- 


-1.4%) 
-1.5%) 



Table 12. Same as Table 11 for the ratio of 14 over 8 TeV cross sections at the LHC. 
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We have checked that the results computed at NNLO+NNLL with NNPDF2.1 are 
fully consistent with the computation of Ref. [84], performed at NLO+NNLL accuracy. 
This proves the reliability of the theoretical uncertainty on the cross section ratios for 
ti production presented in [84], and underscores the stability of this ratio under higher- 
order corrections. Indeed, central values and PDF uncertainties are unchanged, while scale 
uncertainties are further decreased. For instance, at NLO+NNLL the scale uncertainty in 
the 14/8 ratio was about 1%, while with the NNLO+NNLL computation it decreases down 
to 0.3%. 

While no measurement of the ratio between 8 and 7 TeV is available, we can estimate 
the expected precision from the absolute cross-section measurements of Table 2 by assuming 
that all systematic uncertainties (but not the luminosity) are fully correlated between 8 
and 7 TeV. This leads to 

(Atlas) 



4hC8/7(**) = 1-36 + 0.11 pb (8%) 

(CMS) , 
LHC8/7'^ 



(rf^lUtt) = 1.40 + 0.08 pb (6%) . (5.1) 



Note that while ATLAS is higher than CMS at 8 and 7 TeV, this trend is inverted in the 
cross section ratio. 

If no such correlation between the systematic uncertainties is assumed we get 

(Atlas) 



4hc8/7(**') = 1-36 ± 0.20 pb (15%) , 



a 



(CMS) 
LHC8/7 



(tt) = 1.40 + 0.11 pb (8%) , (5.2) 



which illustrates the importance of maximizing the cancellation of systematics. This re- 
quires a dedicated strategy for the ratio measurement, rather than simply combining data 
at two different energies. 

We compare graphically the theoretical predictions for the ratios between 8 and 7 TeV 
and 14 and 8 TeV in Fig. 10. In the first case, we also show the estimate of the ATLAS 
and CMS results, assuming full correlation of the systematic uncertainties. Then we show 
in Fig. 11 the dependence on as{Mz) of the NNLO+NNLL cross-section ratios for each of 
the various PDF sets. The slope with ag for the prediction of each of the PDF sets varies 
more than for the absolute cross sections. Fig. 3. However, with the assumed uncertainty 
on a^, the contribution 5a^ to the total theory uncertainty of the cross-section ratios is 
much smaller than the PDF uncertainty. 



6 Hypothetical fourth-generation heavy quark production at the LHC 

Following Ref. [8], we provide also the total cross section at NNLO+NNLL accuracy for 
a pair of hypothetical heavy fourth-generation quarks, belonging to the fundamental rep- 
resentation of SU{2>). Such new massive fermions arise naturally in BSM theories with 
strongly-coupled dynamics. We denote these hypothetical heavy fourth- generation quarks 
by T. Our aim is to assess the scale and PDF uncertainties affecting the QCD contribution 
to the production of such heavy fermions using the most up-to-date theoretical inputs. We 
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Figure 10. The best predictions from eaclr PDF set for tlic LHC ratios between 8 and 7 TeV (left 
plot) and 14 and 8 TeV (right plot). The error bars correspond to the total theoretical uncertainty. 
For the 8/7 ratio, the ATLAS and CMS experimental results have been obtained taking the ratio of 
the central data and assuming that experimental systcmatics (but not luminosity) is fully correlated 
between different beam energies. 
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Figure 11. The theoretical predictions for the cross section ratios between different LHC beam 
energies for the the various PDF sets, as a function of the strong coupling constant as{Mz)- 



use the MSTW2008 NNLO PDFs as input in the computation, and provide predictions for 
masses niT in the range between 200 GeV and 1.3 TeV. The PDF and scale uncertainties 
are defined as in Sect. 3, with the difference than now the central renormalization and 
factorization scales are set to the heavy quark mass, fip = fJ-R = "niT- 

The numerical results for the total cross sections and associated scale and PDF uncer- 
tainties have been tabulated and they are available in the source of the arXiv submission 
of this paper. In Fig. 12 (left) we show the production cross sections at LHC 7 and 8 TeV 
as a function of rriT, where the uncertainty band is the linear sum of scale and PDF uncer- 
tainties. In Fig. 12 (right) we also show the relative PDF and scale uncertainties, at the 
LHC 8 TeV. We notice that for large heavy fermion masses mT ^ m^, PDF uncertainties 
become the dominant source of theoretical error. 

^The data files are heavyf ermion_7tev . data and heavyf ermion_8tev . data for 7 and 8 TeV results 
respectively. 
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Fourth-Generation Heavy Quark Production at ttie LHC Fourtli-Generation Heavy Quark Production at the LHC 8 TeV 




Mr ( GeV ) Uj ( GeV ) 



Figure 12. Left plot: the cross sections for tlie production of a TT pair of liypotlietical 
heavy fourth-generation quarks at the LHC 7 and 8 TeV, computed at NNLO+NNLL with 
MSTW2008NNLO, as a function of the lieavy quark mass m^. The uncertainty band is the linear 
sum of PDF and scale uncertainties. Right plot: the relative PDF and scale uncertainties as a 
function of mr for LHC 8 TeV. 

7 Summary and outlook 

In this work we have presented a detailed assessment of the theoretical systematic uncer- 
tainties that affect the total top quark pair production cross section at hadron colliders. 
We have also compared the theoretical predictions with the most recent experimental data 
from the Tevatron and the LHC 7 and 8 TeV, and provided predictions for LHC 14 TeV. 

For our analysis we have used the most precise perturbative and non-perturbative 
information available, in particular the recently computed NNLO+NNLL cross sections and 
the most up-to-date NNLO PDF sets. Thanks to the significantly reduced scale dependence 
of the full NNLO calculation, now the three main sources of theoretical uncertainty, scales, 
PDFs and the top quark mass, are roughly of the same size, around 2-3%. We find that 
most PDF sets agree both in their central values and in the size of their PDF uncertainties. 

Given the high accuracy of the perturbative tt cross section, we have demonstrated that 
the available data provide a strong constraint on existing NNLO PDF sets, in particular 
on the large-x gluon PDF. The inclusion of Tevatron and LHC top quark data in the 
NNPDF2.3 set results in a reduction on the large-x gluon PDF uncertainty, which has 
further important phenomenological consequences. We have explicitly studied two such 
cases, both relevant for BSM searches: high mass graviton production in warped extra 
dimensions scenarios, and the high mass tail of the top pair invariant mass distribution. 
We have also given updated predictions for the pair production of a hypothetical fourth- 
generation heavy quark T, which should be helpful in interpreting the results of ongoing 
experimental searches [85-87]. 

For completeness, we have also provided up-to-date predictions for the cross-section 
ratios of at different LHC center-of-mass energies. 

To conclude, we emphasize that the availability of the NNLO calculation for the total 
tt cross section makes top pair production the first hadron collider process which is both 
directly sensitive to the gluon PDF and can be included consistently into a NNLO global 
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PDF analysis without any approximation. This leads us to believe that top quark data 
will be an important ingredient of future PDF fits, especially as more precise data on the 
total cross section and on various differential distributions becomes available. 
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